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Abstract 

We calculate the topological charge density of SU(N) lattice gauge fields for values of 
N up to N = 8. Our T ~ topological susceptibility appears to approach a finite non- 
zero limit at N = oo that is consistent with earlier extrapolations from smaller values of 
N. Near the deconfining temperature T c , we are able to investigate separately the confined 
and deconfined phases, since the transition is quite strongly first order. We find that the 
topological susceptibility of the confined phase is always very similar to that at T = 0. By 
contrast, in the deconfined vacuum at larger N there are no topological fluctuations except 
for rare, isolated and small instantons. This shows that as iV — > oo the large-T suppression of 
large instantons and the large- N suppression of small instantons overlap, even at T ~ T c , so 
as to suppress all topological fluctuations in the deconfined phase. In the confined phase by 
contrast, the size distribution is much the same at all T, becoming more peaked as N grows, 
suggesting that D(p) oc S(p — p c ) at iV = oo, with p c ~ l/T c . 



1 Introduction 



The finite-temperature deconfinement transition of SU(N) gauge theories has, recently, been 
studied numerically for gauge groups ranging from N = 2 to N = 8 (HE]- The transition is 
first order for N > 3 and, as with much of the other physics of these gauge theories [3], the 
approach to the N = oo limit appears to be rapid. 

As part of this ongoing study [4] we have investigated the topological properties of the 
gauge fields in the neighbourhood of the transition. Since for iV > 4 the transition is quite 
strongly first order, it is possible to do this for the confined and deconfined phases separately, 
at exactly the same temperature, when T ~ T c . What happens as iV — > oo is of particular 
interest and what we find is that the topological properties in the confined phase are the same 
for all T as at T = 0, while in the deconfined phase there are no topological fluctuations at 
all, even at T ~ T c . Since at N = oo QCD (with m q > 0) and quenched QCD are the same, 
this conclusion also applies to the deconfined phase of QCDtv=oo- Moreover, as we shall see, 
the suppression of topological fluctuations in the deconfined phase, and the corresponding 
restoration of the t^i(l) symmetry, is very rapid in N (probably exponential). 

We have also performed some SU (6) and SU(8) calculations at T = in order to obtain 
the string tension, a, which then provides a scale in which to express the calculated values of 
the deconfining temperature T c . We have, at the same time, calculated the topological charge 
and this allows us to extend earlier T = calculations, both of the topological susceptibility 
[3J IS] and of the size distribution of the topological charges E|. This size distribution 
becomes narrower as N increases suggesting that it might be tending to a simple ^-function. 

2 Some expectations at large N 

We review here some expectations that one has about topology in the large- N limit, that 
follow from general counting arguments jjj. 

2.1 The topological susceptibility 

The topological susceptibility is defined as \t = (Q 2 )/V, where V is the volume of space-time 
and Q is the total topological charge. Since fluctuations are suppressed as iV — ► oo, so that 
one can imagine N = oo physics being given by a single (gauge-invariant) Master Field jHj, 
and since (Q 2 ) is a measure of fluctuations, the leading-order value of xt should vanish and 
one needs to ask at what order in 1/N one expects it to appear. 
We can write 

(Q 2 ) = £<Q(*)Q(v)> = V E<QOr)Q(0)> (1) 

x,y x 

and we expect (Q(x)Q(y)) to factorise up to relative corrections of 0(1/N 2 ) 

(Q(x)Q(O)) N =°° (Q(x))(Q(0)) + f(x)0(N<- 2 ) = (Q(0)) 2 + f(x)0(N<- 2 ) (2) 
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where we take the ^-dependence of the leading factorised term to be N^, and we use the fact 
that (Q(x)) = (Q(0)) by translation invariance. Putting this together and using (Q(0)} = 
(Q) /V we obtain 

^ = e{^} 2 + E mo(Nn = ^ + e mo(N^). ( 3 ) 

^ x y " ) x v x 

In order to determine ( let us suppose for the moment that we have added a #-term to the 
action, i.e. 5S = i6Q, and that we are working with a generic nonzero value of 9 so that 
(Q) ^ 0. We have 

{Q) = ^hLZ{0)=i±e(6)V (4) 

where we have used Z = exp — eV where e is the vacuum energy per unit volume. Now we 
expect from the usual large- N counting arguments JH] that e oc iV 2 and that a smooth large- N 
limit is reached if one keeps 9/N fixed i.e. 

e (0) = N 2 h{6/N). (5) 

Plugging this into eqn(0J) and using the notation tp = 9/N, we immediately see that 

(Q) = = NVi ip h &) ( 6 ) 

i.e. (Q) oc NV. Thus ( = 2 and we see from eqn® that for any ip ^ 

Xt = {iQ ~y Q)?) oc O(V )O(N°) (7) 

where we have generalised the definition of Xt to allow for (Q) ^ 0. In obtaining eqnfJTJ) we 
use the fact that J2x f( x ) = 0(V°) which follows from the fact that the mass of the lightest 
glueball with the quantum numbers of Q(x), i.e. J PC = h , is non-zero. By continuity we 
expect that at 9 = eqn(J7|) will be valid and that Xt oc 0(N°). 

We therefore expect that Xt should have a finite non-zero limit as N — * oo, even though this 
is really an 0(1/N 2 ) correction in that limit. This is of course nothing but the conventional 
expectation, needed to explain the mass of the rj [TQ1 HI] • Lattice calculations [3 E] support 
this expectation and we shall provide further evidence in this paper. 

2.2 The topological charge 

We have noted above that for 9 ^ we have 

(Q) 



V 



oc N. 



This is not unexpected: instantons are SU(2) objects and (in the semiclassical limit) can 
be simultaneously placed in any of the ~ N/2 non-overlapping subgroups of SU(N) without 
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interacting - even if they are centered at the same point in space-time. At 9 = 0, where we 
work, one continues to expect O(N) charges in every finite region of space-time, although 
now they will largely cancel in (Q). Since the general arguments given in Section \2. II imply 
that (Q 2 ) oc iV° rather than oc iV 1 , this cancellation is clearly much stronger than the usual 
'random walk' of a dilute gas. 

There are detailed variational treatments of the instanton partition function that explicitly 
show how this cancellation can happen ITS] . It is instructive to sketch the core of the 
argument. Since the (anti)instantons can be thrown into oc N SU(2) subgroups in the space- 
time volume V, we can expect their numbers to be n + ,n_ oc VN and equal on the average. 
That is to say n + + n_ oc iV (we drop the ubiquitous factor of V from now on) while the 
distribution in Q = n + — n_ is flat for Q ^ ^JN . (It is easy to see this in an explicit dilute 
gas calculation.) For the fields not to interact they have to be in subgroups that are exactly 
orthogonal and this is a set of fields of measure zero. Allowing fluctuations away from these 
particular subgroups involves interactions. Let 7 s ,7 a label the interactions between like and 
unlike charges that are in the same volume of space-time. The simplest instanton calculation 
[T2] tells us that the interaction is repulsive with 7 a s oc 1/N. Thus the extra interaction piece 
in the action will be 

SSim oc ]p s n + (n + - 1) + -7 s n_(n_ - 1) + 7a«_n + (9) 

which we can rewrite as 

SSi nt oc i^^7 s (n+ + n_) 2 + ]^-^--i s (n + - n_) 2 (10) 

defining r = 7 a /7s (and suppressing a non-leading term linear in n + + n_). The second 
term provides a factor in the partition function of oc exp — cQ 2 where c = 0(1) since the 
factor of 1/N from 7 S will combine with the usual 1/g 2 factor to give the 't Hooft coupling 
that is kept constant as N — > oo. Of course, for c > we must assume that r < 1 i.e. 
7a < 7s, since otherwise the favoured configuration would clearly have charges of the same 
sign, so cancellations would be suppressed and we would have Q oc N. If all this is the case 
then, given that the partition function is otherwise flat in Q = n + — n_ for Q <C v~N, this 
factor of exp —cQ 2 will ensure that (Q 2 ) is 0(1) rather than O(N). What happens to the term 
oc (n+ + nJ) 2 in eqn(JTUJ) is more difficult to see, and requires the full apparatus of a variational 
treatment ^2]- But the suppression of (Q 2 ) from O(N) to 0(1) is (in this instanton picture) 
simply the result of the fact that it induces a suppression in the action that is quadratic in 
Q because it is quadratic in the n + ,n_. This might appear to be a severe breakdown of the 
dilute gas approximation, but it occurs within the same region of space-time, and still leaves 
open the possibility that a dilute gas/liquid picture can describe the interactions of the net 
charges that are located in different regions of space-time. 

It would clearly be interesting to confirm the above expectation, and to explore how it 
works in the full non-perturbative vacuum of the gauge theory. Unfortunately our lattice 
calculations were not specifically designed to address this issue. Nonetheless we do calculate 
two potentially relevant quantities: the total action and Q mo d = J2 X \Q( X )\ as a function of 
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cooling. If there are O(N) topological charges in every region of space-time, and if there 
is some scatter in the positions and sizes of these charges, then this might be revealed in a 
component of these two quantities that grows with N. Unfortunately our calculations of Q mo d 
have shown nothing significant (which may be due to the coarseness of our lattice spacing) and 
while the average plaquette after a given number of 'cooling sweeps' (see Section 0} is roughly 
independent of N, so that the normalised action is oc iV 2 (recall the factor (3 = 2N/g 2 = 
2N 2 /(g 2 N) in the lattice action - see Section El - and that we keep g 2 N fixed as iV — > oo), 
this is difficult to interpret on its own since we of course expect the full action to show this 
behaviour irrespective of topological considerations. 

2.3 Small instantons at all T 

For large N, instantons are exponentially suppressed in N [T4*j : 

O 8t 2 8tt 2 at 

D(p) oc e~ bl oc e = e (11) 

where Sj is the instanton action, p is the instanton size and A = g 2 N is the 't Hooft coupling 
which is kept constant as N — ► oo. This simple semiclassical argument is incomplete (see e.g. 
[T^J El UH] ) but is essentially correct for topological charges that are small enough and using 
the one- loop expression for g 2 (p) we expect that 

D(p) oc e-* P «° p 1 ^- 5 . (12) 

Thus we expect that there is some critical size p c such that 

D(p) ; Vp < p c . (13) 

There is some evidence for this from earlier lattice calculations [3]. 

This has an important practical consequence for lattice calculations. Since the lattice 
Monte Carlo is a local process, in that the elementary step consists of changing one link 
matrix at a time, the value of Q can only change in a sequence of Monte Carlo generated 
gauge fields if an instanton gradually shrinks until it disappears within a lattice hypercube 
(or the converse process). For small lattice spacings this requires the presence of lattice fields 
containing an instanton of size a < /5 < p c . From eqn()13|) the probability of such a lattice 
field — > exponentially fast in N as iV — > oo. Thus at fixed a topology suffers a critical 
slowing down in N, i.e. the value of Q ceases to change and we cannot calculate (Q 2 ) in the 
usual fashion. (Of course, fluctuations corresponding to topological charges of opposite sign 
can still appear and disappear and in principle one could calculate (Q 2 ) from the charge Q 
within a large sub-volume of a very much larger space-time volume - with the latter chosen 
large enough that the constraint of a fixed total topological charge is negligible.) In practice 
this means that in our SU (8) calculations, the value of Q ceases to change as we reduce a and 
we are unable to perform a continuum extrapolation for the susceptibility. Instead we will 
perform a comparison of the susceptibility at a fixed not-too-small value of a. 
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2.4 Large instantons at high T 



At high T, chromo-electric fields are screened on a scale characterised by the electric screening 
mass m e \. At high enough T the characteristic expansion parameter, g 2 {T), is small and it 
suffices to calculate m e \ to leading order in perturbation theory, whereupon one finds flj 

™ 2 el (T) = j 9 2 (T)T 2 . (14) 

Here the power of T follows on simple dimensional grounds, and the coupling appears in the 
familiar combination, g 2 N, which is kept constant as iV — > oo. The chromoelectric fields in 
the core of an instanton are coherent over its size p. On the other hand they will necessarily 
be screened over a distance 0(1 /m e i). This means that instantons of size p 3> l/m e i 1/T 
will be strongly suppressed, and we expect the suppression to be a function of pT. A careful 
calculation [H| shows that the instanton density, D(p), is modified as follows for npT 3> 1 
and g 2 (T) < 1: 

2N 

D(p, T) = D(p, 0) exp(-— {vrpT} 2 - 7 (pT)) (15) 

where j(pT) is a weak function of its argument, which we can neglect for our purposes. We 
note that the dominant suppression in eqn ()15J) is a function of pT as expected. It also contains 
a factor of N, which one would expect for the same reason that it appears in the classical 
action in eqn lfTTj) . 

When valid, eqn (fT5|) implies that topological fluctuations on size scales p ^> 1/nVNT 
will be heavily suppressed, and that at fixed p the suppression will be exponential in N. So 
what is its domain of validity? The derivation of electric screening requires that we be in 
the deconfined phase and the leading-order perturbative calculation is only guaranteed to be 
valid for T ^> T c . So an interesting question is down to what value of T does this calculation 
hold? The most extreme and yet most elegant scenario is one in which it is valid at all T 
in the deconfined phase and that as N — > oo the small instanton suppression described in 
the previous section, overlaps with this larger instanton suppression, so that all topological 
fluctuations vanish exponentially with N at any value of T, in the deconfined phase. This, as 
we shall see below, is indeed what appears to happen. 



3 The lattice calculation 

The lattice calculations have been described in detail elsewhere PjJ |2] and the calculation of 
the topological charge density is precisely as in [3]. We shall therefore restrict ourselves to a 
minimal description here and refer the interested reader to those earlier papers. 

We work on periodic hypercubic lattices . We label the lattice spacing by a and the size, 
in lattice units, by L\L t . For L s sufficiently large this corresponds to the field theory at a 
finite temperature 
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If L s ,L t ^> l/aT c we refer to the system as being at T = 0, although this is of course not 
exactly the case. The variables are SU(N) matrices assigned to the links of the lattice. We 
employ the usual plaquette action jHj which contains the factor 

2N 

P = — (17) 

g 2 

Here g 2 is the bare coupling defined at the cut-off scale a, and by varying the value of (3 we 
vary the value of a. If we tune a — > we recover the continuum field theory in Euclidean 
space-time, if T = 0, or the statistical physics of the theory at temperature T, more generally. 

We evaluate the lattice Feynman Path Integral/Partition Function using standard Monte 
Carlo techniques. At a fixed value of (3 we can calculate the T = string tension in lattice 
units, a 2 ((3)a, and its value allows us to express the lattice spacing a(/3) in physical units. If 
we increase /3 on a lattice with L s ^> L t we will find a deconfining transition at some (3 = (3 C . 
The corresponding deconfining temperature is 

a({3 c )T c = -L (18) 

Dividing this by the T = value of a(f3 = /9 c )a/ct we obtain a value of the dimensionless 
ratio T c j ^Jo which differs from the continuum value by lattice corrections of order a 2 ((3 c ). For 
larger N we have typically calculated (3 C for L t = 5,6,8 and this allows us to make the desired 
a — > continuum extrapolation. 

We calculate the topological charge density using one of the standard methods 01111 • For 
smooth lattice gauge fields one can define a lattice topological charge density Ql(x) as follows: 

Ql&) = ^- 2 ^u pa ReTr{U, v {x)U pa {x)} a -=% a A Q(x) + 0(a e ) (19) 

where Q(x) is the continuum topological charge density. To obtain a smooth lattice field from 
a rough Monte Carlo generated lattice field we apply an iterative 'cooling' procedure [HI EE]. 
Using eqn(Hnj) on the cooled lattice field we obtain its total topological charge Ql = J2 X Ql{x). 
This will not be an integer because of the lattice corrections in eqn (H9|) . but can easily be 
rounded to the appropriate integer. We can also attempt to identify individual instantons by 
looking for peaks in Ql(x) and applying the classical continuum relation 

6 

Qpeak 5 J - (^0) 

to extract the instanton size, p, from the peak height. In this way we can calculate an instanton 
size density, D(p). 

If we use a modest number of cooling sweeps (typically 10 or 20) the calculated topolog- 
ical charge is accurate except in the presence of very narrow charges whose interpretation is 
naturally ambiguous on a lattice. One can see this by comparing with a fermionic calculation 
of Q using fermions that satisfy the index theorem . The topological charge density, on the 
other hand, is steadily distorted by cooling and so the information one obtains on D(p) is less 
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reliable. (Although here too there is some evidence from using massless fermions as a probe, 
that the distortions are not large for a small number of cooling sweeps |19j .) Together with 
the crudeness of eqn ()20)l as a pattern recognition algorithm for a dense 'instanton' ensemble 
(see e.g. j2H] for more sophistication) it is clear that our estimate of D(p) must be regarded as 
semi-quantitative. In particular, while small instantons are easy to identify, since the peak is 
large, very large instantons correspond to very small and smooth peaks which might be rem- 
nant artifacts of the cooling procedure. In this paper we will make some progress in dealing 
with this latter uncertainty. 

4 Results 

4.1 Topological susceptibility at T=0 

In our T = 0' SU(8) calculation, it is only for the coarsest lattice spacings that Q changes 
sufficiently frequently that we obtain a usefully accurate value for (Q 2 ). This means that we 
cannot reliably extrapolate the susceptibility to the continuum limit and compare its value 
with continuum values previously obtained for N < 5 in |3j and for iV < 6 in [H] as we would 
ideally wish to do. Instead we will calculate the dimensionless ratio of the susceptibility, a 4 Xt, 
and the string tension, a 2 a, i.e. axJ A ja^/o = Xt^/V* 7 ' a ^ a fi xe d value of a, chosen so that 
we are able to obtain an accurate SU(8) value. We do so for iV = 2,3,4,6,8. The large- iV 
counting for the lattice gauge theory is the same as for the continuum gauge theory [2I 5 , and 
so one can perform an extrapolation to N = oo, using corrections that are powers of 1/N 2 . 

If we want to use the 'same' value of a in different SU(N) gauge theories, we need to 
measure it in units of some physical quantity that has a smooth large- N limit. We shall 
choose to fix the value of a in units of the deconfining temperature, T c . To be precise, the 
value of a that we choose corresponds to a = 1/5T C . That is to say, we perform our calculations 
at the value of (3 at which the gauge theory on an L 3 5 lattice (with L ^> 5) will undergo a 
deconfining transition. We could have used some other physical quantity, such as the string 
tension, to fix a and the resulting values of (3 would have differed, because T c /y/cr does vary 
with N. But because such physical ratios possess a smooth large- N limit, any difference will 
correspond to an change in the coefficients of terms that are higher order in 1/N 2 . 

In Table n w e list our results for (Q 2 ) for iV = 2,..., 8. We show the lattice sizes, the 
string tensions, and the values of f3 for each value of N. We also show for each N the critical 
value of /3 at a lattice spacing a = 1/5T C . We see that the value of (3 we use is always, within 
errors, equal to (3 c (a = 1/5T C ), i.e. the lattice spacing is constant across all N when expressed 
in units of T c . At a finite lattice spacing different methods of calculating Q will differ. The 
value we show is obtained after 20 cooling sweeps. Q r is the real number obtained directly 
from summing our twisted plaquette operator, while Qi has been rounded to the appropriate 
neighbouring integer value. We show in Table El the resulting values of the dimensionless ratio 
Xt A [ /V a that we wish to compare across N. 

In Fig. ^ we plot the values of Xt^/V* 7 against the expected expansion parameter 1/N 2 . 
To fit the iV-dependence we use as few correction terms as possible. However a simple 1/N 2 
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correction does not provide an acceptable fit, if we insist on fitting all the way down to N = 2 
(or even only down to N — 3). We therefore introduce an additional 1/iV 4 correction and this 
enables us to obtain an acceptable fit. In that case our best fits are 



Xt 



0.397(7) + ^ -i^il using Q, 
0.382(7) + 2^131 usmg Q 7 



(21) 



with a goodness of fit x 2 — 1-0, 1.3 per degree of freedom respectively. (The difference between 
the two fits is an 0(a 2 ) lattice correction and should disappear in the continuum limit.) As 
we see from Fig. the convergence to the N = oo limit is rapid, and our values are accurate 
enough, and extend over a large enough range of N, to make this a reasonably compelling 
conclusion. 

As an aside we note that the coefficient of the l/N 4 term in eqn (f2~T|) is quite large compared 
to the other coefficients, and one might be tempted to ascribe this to the fact that at this 
lattice spacing we are very close to the bulk transition 013111. However, although the original 
a = large- N calculations of Xt^/V a were able to interpolate within errors down to AT = 2 
with just a simple 0(1/N 2 ) correction, and therefore the coefficient of the 0(1/N 4 ) term could 
not be usefully estimated, there are now much more accurate a = calculations [S] from which 
we can extract the coefficient of the l/N 4 term in the continuum limit and we find that it has 
a similar size to the one we have found at a = 1/5T C , although with the opposite sign. 



4.2 Topological fluctuations at T ~ T c 

When the deconfining transition is first order, as it is for AT > 3, and when T is sufficiently 
close to T c , the system will repeatedly tunnel between confined and deconfined phases. At 
T ~ T c the tunneling probability — > as V — > oo, and what that means for our Monte Carlo 
calculation is that for large V we have very long sequences of lattice gauge fields that remain 
in a single phase. Thus we can calculate the topological charge separately in the confined and 
deconfined phases when T is near T c . Indeed there exist calculations of this kind in SU (3) (see 
[22] for early examples) which show that at T ~ T c the topological susceptibility is markedly 
suppressed in the deconfined phase as compared to the confined phase. In this Section we 
extend these calculations to A^ > 3 so as to learn what happens as A" — > oo. 

In Table |3]we present some results on the value that Xt takes in the confined and deconfined 
phases at T ~ T c , for various values of N. The calculations are at a finite lattice spacing and 
the value of a is the same for all A^ when expressed in units of T c : specifically aT c ~ 5. 
We present the susceptibility in the form of the dimensionless ratio Xt/& 2 , which possesses a 
simple interpretation in the limit of a dilute gas of topological charges: 

Xt = (Q 2 ) = (Ni) (22) 

a 2 {aL s } 3 aL t a 2 {aL s } 3 aL t a 2 

where Nj is the total number of topological charges in the field on the L\L t lattice. Thus the 
quantity Xt/& 2 measures the total number of topological charges per unit (four)volume, where 
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the volume is measured in units of the confining string tension. (Recall that in QCDtv =3 , one 
has yja ~ 440MeV ~ 0.45fm.) This interpretation is of course only exact in the dilute gas 
limit, which is certainly not valid in the confining phase of the SU(N) gauge theory, although, 
as we shall see below, it does appear to be the case in the deconfined phase. 

Since we are at a finite a, different ways of calculating Q will differ - typically by 0(a 2 ) 
lattice artifacts. The results we present in the remainder of this sub-section are based on the 
real- valued topological charge, Q r , calculated after 20 cooling sweeps. The results one would 
obtain using the integer value topological charge and/or fewer, say 10, cooling sweeps lead to 
the same conclusions. 

For SU (8) we generated two sequences of 50,000 fields on 12 3 5 lattices with one starting in 
the confined phase and one in the deconfined phase, and with the inverse coupling (3 chosen so 
that T = l/5a(j3) ~ T c . The spatial volume is large enough that not only are there no phase 
transitions but there is no sign of even a partial tunneling. (We use the time-like Polyakov 
loop averaged over the lattice as an order parameter to discriminate phases.) Thus there is 
no ambiguity in calculating quantities separately in the confined and deconfined phases. For 
SU(6) we do the same (with sequences of 100,000 fields) on 16 3 5 lattices, for S"?7(4) (with 
sequences of 25,000 fields) on 32 3 5 lattices and for SU(3) on 32 3 5 lattices. In the last case the 
volume is small enough that we observe several tunnelings between the phases, but these are 
sufficiently well-defined that there is no significant ambiguity in separating the phases from 
each other. That one can work on smaller volumes as N f is mainly due to the fact that the 
surface tension of the bubble separating the confined and deconfined vacua increases with iV 
[2] so that for a given volume V the tunneling probability — > as N — > oo at T = T c , and so 
we have longer sequences of fields in which the system remains in a single phase. Finally we 
remark that while we have not included SU(2) here, since its transition is second order, one 
also finds in that case a suppression of Xt for T > T c 

For comparison we also show in Table|3]the 'T=0' value of Xt/& 2 , as obtained from Tabled 
(or a slight extrapolation thereof). Comparing the various values in that Table, we see that 
Xt/cr 2 remains approximately constant in the confining phase at all T, but that it is strongly 
suppressed in the deconfined phase. Moreover the suppression becomes rapidly more severe 
as N grows. To illustrate this we show in Fig. |2] how the ratio of Xt i n the deconfined phase 
to its value in the confined phase varies with N at T = T c . It seems quite clear that in the 
deconfined phase Xt(T = T c ) — ► as N — > oo. 

In Tables 0] to [7|we list some further calculations, some at higher T. (The values marked 
by a * in Tables H3 to [7| have been obtained from sequences of fields containing sufficient 
tunnelings for there to be some danger that the quoted error is underestimated.) Together 
with our other calculations, these provide more information on the T dependence of Xt- As 
an example, we plot in Fig. 01 the T dependence of Xt/& 2 for SU(6). We see that the value 
of Xt in the deconfined phase shows a strong temperature dependence, vanishing rapidly as T 
grows. 
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We conclude that at iV = oo a rather simple picture emerges, 




const 0) 







deconfined VT, 



confined VT, 



(23) 



with the suppression in the deconfined phase being very rapid - probably exponential in a 
power of N. This shows that the the rapid instanton suppression, for which we have reliable 
semiclassical arguments only at high-T where g 2 {T) is small, in fact occurs at all values of T, 
even down to T = T c , and indeed for all values of the instanton size p. The confined phase, 
on the other hand, appears to be impervious to the large-p suppression, even at T ~ T c . To 
make this interpretation convincing, we need to estimate the instanton size distribution, and 
this we shall now do. 

4.3 Instanton size distributions 

To obtain information on the instanton size density, D(p), we smoothen ('cool') the lattice 
fields, so that our lattice topological charge density operator should approximate the con- 
tinuum one, and we then locate the peaks in that density, Q(x). We identify each of these 
with an instanton, whose size p is given by eqn (J2(J|) . (As described in Sectional) This pro- 
cedure is clearly unambiguous for small instantons, which correspond to large peaks in Q(x), 
but becomes increasingly more ambiguous for larger instantons that correspond to very low 
and smooth bumps in Q(x). To help avoid mere fluctuations in the instanton profile being 
misinterpreted as extra instantons, we do not count peaks that are within a distance of 2a 
of a larger peak. This algorithm is of course very crude (as emphasised in Section EJ) but 
unfortunately we are not aware of a really reliable alternative. We shall therefore focus on the 
qualitative features of our calculated size densities. 

4.3.1 T=0 

The T = size density was calculated for 2 < N < 5 in j3J. Here we shall extend the 
comparison to iV = 8. However none of our lattice spacings is as fine as that used in Fig. 13 
of [3] and this means that our comparison will be more affected by lattice spacing artifacts. 

We start with the calculations listed in Table [TJ In each case we calculated the topological 
charge density after 20 cooling sweeps and the resulting number densities are shown in Fig. |3J 
What we plot is the number of charges per 10 4 lattice fields with a = 1/5T C . (If one wishes 
one can translate the volume into more 'physical' units using the listed values of the string 
tension, a^/a, and using the real world value of ^/o ~ 440MeV.) The histograms are for 
numbers of instantons in bins of p that are of size Sp = 0.25a. Each point is placed at the 
average value of p within that bin, rather than in the middle of the bin, which means that the 
points for different N are slightly shifted (in p) with respect to each other. 

The value of a is very coarse and this means that the effective cut-off induced by the 
cooling, which here we can estimate to be at about p ~ 2 — 2.5a from the 577(2) density, is 
uncomfortably close to the typical instanton size. Nonetheless some qualitative features are 
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clear. First is the steepening of the distribution at small p. This is expected theoretically and 
was already observed up to N = 5 in [3] . More interesting, and only becoming apparent from 
the greater range of N available to us in the present calculation, is that the whole distribution, 
at both small and large p, is becoming much more peaked as N grows. Indeed, if we take the 
distributions in Fig. H] seriously, we infer that 

lim D(p) = c(N)5(p - Pc ) ; p c ~ 1 (24) 

The main uncertainty is to do with the reliability of our identification of the larger topological 
charges. However since the size distribution narrows with increasing N, the relevant range 
of p comes closer to the average, and so our calculation should be more reliable at larger N. 
Some explicit checks showing that our conclusions are robust against the detailed cuts of our 
algorithm would be useful, but have not been performed in the present calculation. Equally, 
a systematic finite volume study would be useful. For the case of SU(3) there exist such finite 
volume studies (e.g. [2U]) and for higher iV our present study has larger spatial volumes at 
T ~ T c , which suggest that finite volume corrections are unimportant. 

The integrated number of topological charges also appears to increase with N, as shown 
in Fig. 03 One might imagine that one is seeing the theoretically expected growth that was 
discussed in Section 12.21 However one needs to be cautious. The volumes at different iV that 
are used here are only constant if expressed in units of T c . If instead they are expressed in 
units of, say, ^/a most of the variation that we see in Fig. |S]for N < 4 disappears. This is 
because the volume goes as the 4'th power of the scale. It is only for iV > 6 that we are close 
enough to iV = oo for this ambiguity to become negligible. 

A study similar to that in Fig. 0] but for a smaller value of a would be very useful. Unfor- 
tunately we have only performed such a calculation for SU(8). In Fig.Elwe show the resulting 
SU(8) instanton size distribution obtained on a 16 4 lattice with a ~ 1/8T C . To compare with 
distributions at smaller N we use data on 16 4 lattices obtained in These are mostly at 
rather different values of a and we rescale the sizes and volumes to a ~ 1/8T C in each case. 
The resulting distributions are plotted in Fig. |U1 The cooling/lattice cut-offs are now far from 
the typical instanton size, and we have an extended region of p in which we can observe the 
dramatic and expected suppression of small instantons. Again we see a rapid narrowing of 
the distribution as N grows, consistent with eqn (J24j) . The total number of charges appears to 
grow with N, but once again there are various reasons for considering this to be a much less 
robust observation. 

4.3.2 T > T c 

We saw in Section 14.21 how, in the deconfined phase, the fluctuations of Q are rapidly sup- 
pressed as N grows. In this section we shall look at the size distribution, D(p), in that phase. 
In doing so we shall continue to use the relation eqn (|2U|) to extract p, although for larger in- 
stantons at high T the solution changes (becoming periodic), as described in J7j. The reason 
is that both in theory and in practice such periodic instantons are highly suppressed at larger 
N. 
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Before looking at the deconfined phase we briefly look at the instanton density in the 
confined phase, and compare densities at T ~ and T = T c . We focus on SU(8). The 
T = T c calculation has been performed on a 12 3 5 lattice at (3 = 43.965 and is compared in 
Fig. to the T ~ calculation on a 10 4 lattice at /3 = 44.0. In the Fi srure sizes have been 
rescaled so that the numbers are per 10 4 volume at (3 = 44.0. What we see is that there is no 
significant difference between the densities at T ~ and at T ~ T c , just as we saw in Table El 
no significant difference in the susceptibilities. The same pattern appears at other values of 
N. We conclude that the topological content of the confined phase is very much the same at 
any value of T. 

Turning now to the deconfined phase of the SU(8) theory, we plot in Fig. El the size 
distribution obtained on 500 12 3 5 lattice fields that are in the deconfined phase and generated 
at (3 — 43.965, i.e. at T ~ T c . We see that the distribution is entirely located at very large 
values of p (compare Fig. [7| which is in the confined phase) except for a very few peaks (4 out 
of about 6000) that lie in the range p G [2, 3.5]. 

To interpret all this we take these 500 fields and use only the largest positive and negative 
peaks in each field. We then take separately the 4 configurations with Q ^ and the 496 lattice 
fields with Q = 0, and we look in Fig. El at the size distributions of the two sets separately. 
Moreover in the former case we show separately the size distribution from peaks with the same 
sign as Q and from those with the opposite sign. We see that the latter have large values 
of p, just like the Q = fields. The former, by contrast, provide the four values e [2, 3.5] 
that we noted in the previous paragraph. The interpretation is clear: real topological charges 
are narrow while the other very broad charges are nothing but artifacts of the partial cooling 
and do not correspond to real topological charges at all. Given that the high-T suppression 
we expect theoretically, oc exp(— ^N{npT} 2 ), is more severe at larger p, the observed gap in 
the size distribution makes the interpretation of the very low, broad peaks as artifacts quite 
compelling; as does the fact that these 'large charges' do not contribute anything to the net 
topological charge Q. 

It is interesting to note that for the confined phase on the same lattice at the same T we 
find 7 charges G [2, 3] (with none narrower). This represents the tail of the highly suppressed 
small instantons. What we see in the deconfined phase is consistent with being the same 
tail. In contrast to the confining phase however, where the suppression rapidly relaxes as p 
increases, here there is a cut-off immediately above p ~ 3.5. Thus it would appear that the 
the high-T exponential suppression is already effective down to small sizes p ~ 3.5 and down 
to modest temperatures T ~ T c . 

Our results at other iV show the same features, although the interpretation is simplest in 
the case of SU(8) because here the deconfined configurations with Q ^ are so rare, so that 
the Q = fields are highly unlikely to contain any instanton-antiinstanton pairs that would 
complicate the argument and the size distributions. To illustrate this we show in Fig. El the 
SU(6) deconfined instanton size density. It show the same features as we saw in Fig. |H] except 
that the low-p peak is much larger - as we would expect because the small-p suppression is 
more severe for SU(8) than for SU(6). For comparison we show the confined density at T c 
in Fig. ^2 The total number of peaks for p < 3.0 is roughly the same in both cases showing 
that what is involved is the suppression of charges with larger p. (From the ratios of the 
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densities in the two Figures one can map out the suppression as a function of p.) In Fig. ED 
we show the SU(Q) version of Fig. El Unlike the latter we now have a small low-p peak from 
the Q = lattice fields. This is because, as we can see, the number of lattice fields with 
Q 7^ is much larger and so there is now a finite probability to have lattice fields with well- 
separated instanton-antiinstanton pairs. For comparison, we show in Fig. ED the same plot 
for the confined phase. 

We conclude therefore that in the deconfined phase the severe small-p suppression is 
matched by an even more severe larger p suppression which, already for N = 8 and even 
at T ~ T c , extends down to very small p. It seems clear that as N —>■ oo the deconfined phase 
will lose its topological charges exponentially fast at any given p; and that this will be so for 
all T in this phase. By contrast the confined phase appears to show no variation with T, even 
at T ~ T c . Clearly there the semiclassical prediction is obstructed by the non-perturbative 
fluctuations. We thus have a rather simple picture of the topological properties of the theory 
as N — > oo. 

5 Conclusions 

As we outlined above, there are various theoretical arguments for the large-iV behaviour of 
topological fluctuations in gauge theories, which it is interesting to test in lattice calculations. 
One expects the topological susceptibility Xt to have a finite and non-zero iV — > oo limit (see 
Section EH) , despite the fact that simple arguments (see Section |2~2J| suggest that the number 
of topological charges grows oc N. One expects smaller instantons, up to some critical size, 
to be completely suppressed as iV — > oo (see Section 12. 3|) and, since this is a short distance 
argument, it should occur at all T. Once the temperature T is high enough for the perturbative 
evaluation of fluctuations around instantons to be reliable, one also expects instantons to be 
exponentially suppressed in N{TrpT} 2 (see Section E3J). 

Compared to earlier calculations, our range of N is larger, and we explore both the T and 
the N dependence of the topological fluctuations. Since the deconfining phase transition is 
robustly first order for larger N, we are able to investigate the confined and deconfined phases 
separately at the same value of T ~ T c . 

Our results on the T ~ topological susceptibility (see Fig. HJ) show that the topological 
susceptibility (albeit at a fixed non-zero value of the lattice spacing) does indeed have a finite 
and non-zero N — > oo limit, as suggested by earlier work. The size distribution shows the 
expected suppression of small topological charges, not only at T ~ (see Fig. |U and Fig. EJ) 
but also at non-zero T, in both the confined (see Fig.EJ) and the deconfined (see Fig. IHJ) phases. 
Within the confined phase there appears to be little if any variation with T, at any T, in the 
character or density of the topological fluctuations. In addition we find that in this phase 
the size distribution appears to become more peaked as N grows, suggesting that it becomes 
oc S(p — p c ) at N = oo. 

In the deconfined phase we find a strong suppression of all topological fluctuations which 
becomes rapidly more severe with increasing N. (See Fig. [2] and Fig. 0) We infer that the 
large-T exponential suppression occurs at all T > T c at large N, and that it overlaps with the 
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small-p suppression so that all topological fluctuations are suppressed. 

The low, broad peaks in the deconfined phase (see Fig. EJ) are remnant artifacts of the 
incomplete cooling. They provide us, for the first time, with a measure of the 'background' 
that such artifacts will contribute to all our other size distributions. Of course, these artifacts 
will depend to some extent on the content of the fields, which is not the same for the confined 
and deconfined phases. Nonetheless the fact that these artifacts are at such large values of p 
reassures us that the part of the size distributions that is important to our conclusions about 
the confined phase (in particular the peaking) is robust. 

Our calculations thus support a very simple picture of topological fluctuations in the 
N = oo limit. In the deconfined phase these disappear completely for all T and for all p. In 
the confined phase all fluctuations with p less than some p c ~ 1/T C disappear, and there is 
evidence that the same is true for all p > p c , so that the size distribution becomes a simple 
5-function. 
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N 


P 


p c (aT c = 0.2) 


lattice 


a^Ja 


(Qi) 




2 


2.3715 


2.3714(6) 


12 4 


0.2879(13) 


3.69(13) 


3.39(12) 


3 


5.8000 


5.8000(5) 


10 4 


0.3133(13) 


2.917(79) 


2.469(74) 


4 


10.637 


10.637(1) 


10 4 


0.3254(14) 


3.375(90) 


2.881(81) 


6 


24.515 


24.514(3) 


10 4 


0.3385(15) 


3.48(18) 


2.96(17) 


8 


44.000 


43.98(3) 


10 4 


0.3413(13) 


3.20(42) 


2.68(35) 



Table 1: Fluctuation of the topological charge in various SU(N) lattice gauge theories for a 
fixed value of the lattice spacing a = 1/5T C . Q r is the raw real- valued lattice charge and Qi 
is our best estimate of the integer that corresponds to it. Also shown is the string tension, a, 
the critical value of /3, and the lattice parameters. 



N 


P 


V 4 / / 

Xi, t N° 


Xr't 1 V a 


2 


2.3715 


0.4013(40) 


0.3928(40) 


3 


5.8000 


0.4171(34) 


0.4001(35) 


4 


10.637 


0.4165(33) 


0.4004(33) 


6 


24.515 


0.4034(56) 


0.3875(58) 


8 


44.000 


0.392(13) 


0.375(13) 



Table 2: The topological suceptibility in units of the string tension for the calculations in 
Tabled 







N 


P 


lattice 




T/T c 


confined 


deconfined 


T = 


3 


5.8000 


32 3 5 


0.3133(13) 


1.00 


0.0270(13)* 


0.0141(13)* 


0.0256(9) 


4 


10.635 


32 3 5 


0.3262(15) 


1.00 


0.0260(19) 


0.0070(7) 


0.0257(9) 


6 


24.515 


16 3 5 


0.3385(15) 


1.00 


0.0199(14) 


0.00171(14) 


0.0226(13) 


8 


43.965 


12 3 5 


0.3462(15) 


0.995 


0.0161(28) 


0.00005(3) 


0.0198(26) 



Table 3: Topological susceptibility in units of the string tension calculated separately in the 
confined and deconfined phases, for the values of N, on the lattice volumes and at the values 
of j3 shown. 



SU(8) 


P 


lattice 


a^J o 


T/T c 


Xr,t/o- 2 :conf 


Xr,t/c" 2 :deconf 


44.00 


8 3 4 


0.3413(13) 


1.26 




0.00000(4) 


43.85 


8 4 


0.3630(21) 





0.0191(24) 





Table 4: Some additional SU (8) calculations of the topological susceptibility in units of the 
string tension. 



16 



SU{6) 




lattice 




T/T c 


Xr,t/o- 2 :conf 


X r t/(T 2 :deconf 


24.500 


16 3 5 


0.3414(20) 


0.99 




0.00219(18) 


24.530 


16 3 5 


0.3356(20) 


1.01 




0.00115(10) 


24.515 


10 3 4 


0.3385(15) 


1.25 




0.000012(5) 



Table 5: Some additional 577(6) calculations of the topological susceptibility in units of the 
string tension. 



577(4) 


P 


lattice 




T/T c 


Xr,i/o- 2 :conf 


Xr,t/o" 2 :deconf 


10.645 


12 3 5 


0.3222(16) 


1.01 




0.0058(3)* 


10.642 


20 3 5 


0.3235(16) 


1.01 




0.0060(3) 


10.637 


20 3 5 


0.3254(14) 


1.00 


0.0246(10)* 


0.0071(4)* 


10.635 


20 3 5 


0.3262(16) 


1.00 


0.0247(11)* 


0.0076(4)* 


10.633 


20 3 5 


0.3270(16) 


1.00 


0.0257(8)* 


0.0074(5)* 


10.637 


12 3 4 


0.3254(14) 


1.00 




0.000115(22) 



Table 6: Some additional 577(4) calculations of the topological susceptibility in units of the 
string tension. 



5f7(3) 


P 


lattice 




T/T c 


Xr,i/o" 2 :conf 


Xr,t/^ 2 :deconf 


5.7950 


20 3 5 


0.3163(20) 


0.99 


0.0243(10) 




5.7975 


20 3 5 


0.3148(16) 


1.00 


0.0253(8)* 


0.0144(10)* 


5.8000 


20 3 5 


0.3133(13) 


1.00 




0.0132(5) 


5.8025 


20 3 5 


0.3118(16) 


1.00 




0.0124(5) 


5.8000 


12 3 4 


0.3133(13) 


1.25 




0.00042(4) 



Table 7: Some additional SU (3) calculations of the topological susceptibility in units of the 
string tension. 
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Figure 1: The topological susceptibility, Xt, in units of the string tension, a, plotted versus 
1/N 2 . The calculations are performed at a fixed lattice spacing a ~ 1/5T C . The line is a 
large- N extrapolation that includes 0(1 /N 2 ) and 0(l/iV 4 ) corrections. 
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Figure 2: The ratio of the topological susceptibility, Xt, in the deconfined and confined phases 



at T ~ T c . 
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Figure 3: The topological susceptibility, Xti i n units of the string tension, a, plotted ver- 
sus T/T c for SU(6), plotted separately for the confined, o, and deconfined, •, phases. The 
calculations are performed at a lattice spacing a ~ 1/5T C . 
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Figure 4: The instanton size density, D(p), for N = 2(*), 3(+), 4(o), 6(x), 8(«) on (mostly) 
10 4 lattices with a = 1/5T C . 
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Figure 5: The average number of instantons versus iV at a = 1/5T C on a 10 4 volume, from the 
calculations listed in Tabled 
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Figure 6: The instanton size density, D(p), for N = 2(*), 3(+), 4(o), 8(«) on 16 4 lattices with 
a ~ 1/8T C . 



23 



2.25 - 



D(p) 1.75 



1.5 
1.25 



0.75 
0.5 
0.25 




n 1 1 1 1 1 1 1 r 



•xx 

x 



x 



x 



x 



x 



x 



X. 



X 



X 



J I , 3ft L 



X. 

X 

XX. 

X*v 

J I I W^w 



6 



8 9 10 



P 



Figure 7: The 5*7(8) instanton size density, D(p), at T ~ (•) and at T ~ T c (x) in the 
confined phase, with a ~ 1/5T C in both cases. Normalised to a 10 4 volume. 
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Figure 8: The SU(8) instanton size density, D(p), at T ~ T c in the deconfined phase, with 
a ~ 1/5T C . Normalised to a 10 4 volume. 
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Figure 9: The number of instantons in 500 lattice fields using only the largest positive and 
negative peaks in each field. All in the deconfined phase, at T = T c , on a 12 3 5 lattice in SU (8). 
Separately for the 4 configurations with Q ^ and for the 496 lattice fields with Q = 0. The 
latter (x) have been divided by 100 so as to fit on the same plot. For the former we show 
separately sizes from peaks with the same sign as Q (•) and with the opposite sign (o). 
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Figure 10: The SU(6) instanton size density, D(p), at T ~ T c in the deconfined phase, with 
a ~ 1/5T C . Normalised to a 16 3 5 volume. 
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Figure 11: The SU(6) instanton size density, D(p), at T ~ T c in the confined phase, with 
a ~ 1/5T C . Normalised to a 16 3 5 volume. 
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Figure 12: The number of instantons in 2000 lattice fields using only the largest positive and 
negative peaks in each field. All in the deconfined phase, at T = T c , on a 16 3 5 lattice in 
SU(6). Separately for the lattice fields with Q ^ and for the fields with Q = 0, x. For the 
former we show separately sizes from peaks with the same sign as Q (•) and with the opposite 
sign (o). 
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Figure 13: The number of instantons in 1000 lattice fields using only the largest positive and 
negative peaks in each field. All in the confined phase, at T = T c , on a 16 3 5 lattice in SU(6). 
Separately for the lattice fields with Q ^ and for the fields with Q = 0, x. For the former 
we show separately sizes from peaks with the same sign as Q (•) and with the opposite sign 
(o). 
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